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ABSTRACT 

We evaluate the conductivity tensor for molecular gas at densities ranging from 10 4 
to 10 15 cm -3 for a variety of grain models. The Hall contribution to the conductivity 
has generally been neglected in treatments of the dynamics of molecular gas. We find 
that it is not important if only 0.1 /im grains are considered, but for a Mathis-Rumpl- 
Nordsieck grain-size distribution (with or without PAHs) it becomes important for 
densities between 10 7 and 10 n cm~ 3 . If PAHs are included, this range is reduced to 
10 9 -10 10 cm- 3 . 

The consequences for the magnetic field evolution and dynamics of dense molecular 
gas are profound. To illustrate this, we consider the propagation of Alfven waves under 
these conditions. A linear analysis yields a dispersion relation valid for frequencies 
below the neutral collision frequencies of the charged species. The dispersion relation 
shows that there is a pair of circularly polarised modes with distinct propagation 
speeds and damping rates. We note that the gravitational collapse of dense cloud 
cores may be substantially modified by the Hall term. 

Key words: magnetohydrodynamics - Alfven waves - dust grains - molecular clouds. 



1 INTRODUCTION 

The breakdown of flux-freezing in molecular clouds has gen- 
erally been treated using the concept of "ambipolar diffu- 
sion" (e.g. Spitzer 1978). In this approximation the magnetic 
field is regarded as being frozen into the ionised component 
of the fluid, which drifts as a whole with respect to the 
neutrals. During ambipolar diffusion, the magnetic stresses 
operating on the ionised plasma are transmitted to the dom- 
inant neutral component via collisions. The resulting dissi- 
pation damps MHD waves (Kulsrud & Pearce 1969) and is 
particularly important in shock waves where it provides the 
heating within the shock front (Draine 1980). Ambipolar dif- 
fusion allows the trickling of neutrals through the ions and 
magnetic field towards the centre of a dense core until the 
core becomes gravitationally unstable and collapses (see e.g. 
Mouschovias 1987; Fiedler & Mouschovias 1993). 

The ambipolar-diffusion approximation is valid when 
the magnetic force on each charged particle dominates the 
drag force arising through collisions with the neutrals (i.e. 
the product of the gyrofrequency and the time scale for 
momentum-exchange with the neutrals is much greater than 
unity). This is generally an excellent approximation for ions 
and electrons in molecular clouds, except for densities and 
magnetic field strengths appropriate for protostellar disks 
(Norman & Heyvaerts 1985; Wardle & Konigl 1993). 

Charged grains, however, are partially decoupled from 
the magnetic field by neutral collisions because of their large 



geometrical cross-section. This effect has generally been 
treated by including a separate equation of motion for the 
grains so that their drift speed through the neutrals can be 
calculated (see e.g. Draine 1980). The current density then 
has a component out of the plane containing the magnetic 
field and the electric field in the rest frame of the neutral 
fluid. This component, the Hall current, has generally been 
suppressed in calculations of gravitational condensation via 
ambipolar diffusion (Ciolek & Mouschovias 1993), and in 
shock modelling (e.g. Draine, Roberge & Dalgarno 1983; 
Wardle & Draine 1987; Kaufman & Neufeld 1996). This is a 
poor approximation when the grain drag is important (typ- 
ically at densities in excess of 10 6 cm -3 . In the case of shock 
waves, when this component is retained the shock structure 
is significantly thinner (Pilipp, Hartquist & Havnes 1990) 
and is qualitatively different, exhibiting a twist back and 
forth around the preshock normal (Wardle 1998). 

Early studies assumed that grains were of a single char- 
acteristic size (typically 0.1 /xm). It has since been recognised 
that the MRN grain-size distribution (Mathis, Rumpl & 
Nordsieck 1977) and the presence of polycyclic aromatic hy- 
drocarbons (PAHs) (Leger & Puget 1984) imply that grains 
play a more important role in coupling the magnetic field 
to the neutral gas (Nishi, Nakano & Umebayashi 1991), 
grains and PAHs being the dominant charged species above 
10 r cm -3 . 

Here we examine the conductivity of molecular gas at 
densities relevant to cloud cores, and at higher densities oc- 
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curring during the formation of a protostar and its associ- 
ated disk. We formulate MHD in in terms of a conduc- 
tivity tensor appropriate for a weakly-ionised gas (Cowling 
1957; Norman & Heyvaerts 1985; Nakano & Umebayashi 
1986) and stress that the vector evolution of the magnetic 
field depends on the relative magnitudes of the different 
components. We evaluate the conductivity for single-size, 
MRN, and MRN-plus-PAH grain models in §||, with charged 
particle abundances from Umebayashi & Nakano (1990) and 
Nishi et al (1991). We find that the Hall term is not partic- 
ularly important in the single-size grain model, but in the 
other models becomes important for densities between 10 7 
and 10 11 cm -3 . The dynamical behaviour of molecular gas 
at higher densities is therefore profoundly different from the 
ambipolar diffusion case. By way of illustration we consider 
the propagation of Alfven waves in Slj. The implications of 
our results are further discussed in §fl, and our conclusions 
are summarised in Sfi 



2 FORMULATION 

2.1 MHD in weakly-ionised media 

Our definition of "weakly-ionised" is that the abundances 
of charged species are so low that their inertia and thermal 
pressure is negligible, and that the changes in the neutral 
gas resulting from ionisation and recombination are so small 
that they may also be neglected. 

In this limit, the fluid equations may be written 

dp 

at 

9v , 2 „ J x B 

p— + p{v ■ V)w + c s X7p = — - — 

where c 2 a , the isothermal sound speed, is assumed to be con- 
stant, 



+ V • (pv) = , 



(1) 



(2) 



J = — V x B , 

47T 

dB 



dt 

and 



V X (v X B) - cV X E' , 



V-B = 0. 

Here 



E' = E + v X B/c 



(3) 
(4) 

(5) 
(6) 



is the electric field in the frame comoving with the fluid, 
which is related to the current density by the conductivity 
tensor a: 



J = cr E' . 



(7) 



We shall consider an explicit expression for a in the 
next subsection. For now, we merely note that all of the 
information concerning the charged species is hidden in er, 
and that we are implicitly assuming a prescription for de- 
termining the abundances of the charged species so that the 
conductivity can be calculated. 

At this point we emphasise the role that the conductiv- 
ity plays in determining the evolution of the magnetic field 
through the induction equation (^), which on substituting 
for E' becomes: 



^-=Vx(vxB)-^-Vx[a- 1 .(VxB)]. (8) 

The second term becomes important if the current asso- 
ciated with the gradient in the magnetic field is so large 
that \E'\ is comparable to or greater than \v X B\/c. The 
magnitude of the conductivity determines a characteristic 
length scale (or time scale) below which the field cannot be 
regarded as frozen-in (see e.g. Parker 1979). It is not, how- 
ever, the magnitude on which we wish to focus but rather 
the tensor nature of the conductivity, which affects the di- 
rection of dB/dt. That is, the vector evolution of B from 
an initial configuration depends on the relative magnitude 
of the components of a. 



2.2 The conductivity tensor 

We characterise the various charged species in the neutral 
fluid by particle mass m,j and charge Zje, number density 
rij , and drift velocity through the neutral gas v j , where the 
subscript j denotes different species. Note that we have al- 
ready implicitly assumed charge neutrality, i.e. 



0. 



0) 



Assuming that the fluid evolves on a time scale that is 
long compared to the collision time scale of any type of 
charged particle with the neutrals, each charged particle 
drifts through the neutrals at a rate and direction deter- 
mined by the instantaneous Lorentz force on the particle: 



Z je (E' + ^X B)- lj m 3 v j = 0, 



(10) 



where the third term represents the drag force contributed 
by collisions with the neutrals, and 



7i 



< av >j 
rrij + m 



(11) 



where < av >j is the rate coefficient for momentum transfer 
by collisions with the neutrals and m is the mean neutral 
particle mass. The Hall parameter for species j, 

ZieB 1 



0j = — . 

mjC jjp 



(12) 



determines the relative importance of the Lorentz and drag 
forces in balancing the electric force. 

Eq ([l(^) can be inverted for Vj, which can then be used 
to form an expression for J = e ^ UjZjVj-. 



J = o-nEi, + a ± B X E' ± + a 2 E' ± 



(13) 



where E'^ and E' ± are the decomposition of E' into vec- 
tors parallel and perpendicular to B respectively. The com- 
ponents of <r are the conductivity parallel to the magnetic 
field, 



6C V — > 
CT || = -g2-^ n 3 Z 3pj > 

i 

the Hall conductivity, 



b ^ i + m 



(14) 



(15) 



and the Pedersen conductivity 
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Figure 1. Upper panel: the abundances of charged species when 
grains have size 0.1 fivci. The curves labelled i, e, G + and G~ re- 
fer to ions, electrons and positively and negatively charged grains 
respectively (adapted from Umebayashi & Nakano 1990). Lower 
panel: the field-parallel (cry), Hall (<ri) and Pedersen ((72) com- 
ponents of the conductivity tensor. The ambipolar-diffusion ap- 
proximation is valid for |<ri| «U2< ani. The spike in the curve 
for |iti| at nn ~ 10 10 cm -3 occurs because <ti is negative at lower 
densities and positive at higher densities. 
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°*=b2^TTW (16) 

3 

(Cowling 1957; Norman & Heyvaerts 1985; Nakano & Ume- 
bayashi 1986). We shall find it useful to refer to the total 
conductivity perpendicular to the field, 

<T X = y/ (17) 

Three conductivity regimes are delineated by the mag- 
nitude of the typical Hall parameter, /?, of the charged 
species: 

(i) \j3\ ^> 1, which implies that cm 3> 02 3> and gives 
rise to the ambipolar diffusion regime; 

(ii) <C 1, which yields o\\ w a% ^ |ci|, implying that 
the conductivity is scalar (J ~ o~^E') - the resistive regime; 
and 

(iii) |/3| « 1, the Hall regime. 

The latter has been largely neglected, but we shall see that 
the grain-size distribution implies that it is relevant over a 
broad range of conditions in dense molecular gas. 



Figure 2. The contributions made by different species to the 
components of the conductivity tensor plotted in Fig. |lj Note that 
the positive and negative species contribute positive and negative 
terms to <ri respectively, but that all species contribute positively 
to try and (72. 

3 CONDUCTIVITY IN DENSE MOLECULAR 
GAS 

The conductivity tensor depends on the abundances of 
charged species, which in turn depend on the choice of grain 
model as recombination generally occurs on the surface of 
grains (Nishi et al 1991). We adopt three grain models: 

(i) Single-size grains, characterised by a grain radius 
0.1 /jm and a total grain mass that is one percent of the 
total mass in hydrogen. 

(ii) An MRN model with a power-law distribution of grain 
sizes between 50 and 2500 A: 

= Anna- 3 ' 5 , (18) 

da 

where n(a) is the number density of grains with radii smaller 
than a, and A = 1.5 x 10~ 25 cm" 2 ' 5 (Draine & Lee 1984). 

(iii) An MRN model with an additional component of 
very small grains with a ~ 3 A, representing PAHs with 
an abundance n g /nn = 2 x 10 -7 . 
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To evaluate the conductivity, we obtain the abundances 
of charged species from Fig. 2 of Umebayashi & Nakano 
(1990) and Figs 1 and 4 of Nishi et al (1991). The species 
we include are molecular and metal ions i, electrons e, MRN 
grains (G + , G - ) and PAHs (g + , g~). The rate coefficients 
for momentum transfer by elastic scattering of species j with 
neutrals, < av >j are given by 



< av >i=< av > 9 = 1.6 x 10 9 cm 3 s 1 



„ 10 _i5 2 /128kT c 
< av > e = 1 x 10 15 cm 2 ' 



V 97rm c 



1/2 



and 



< av >o= na 



/128fcT\ 1/2 



(19) 
(20) 

(21) 



V 9irm I 

where T and T e are the neutral and electron temperatures 
respectively. The rate coefficient for the PAHs has been set 
equal to that for the ions as the induced-polarisation cross- 
section dominates the geometric cross-section for charged 
particles with sizes below 10 A. The expressions for < av >c 
and < av > e are valid as long as the grain and electron drift 
speeds are less than the sound speeds in the neutrals and 
electrons respectively. 

To calculate the conductivities, we adopt T — T e = 
30 K, a cosmic-ray ionisation rate per hydrogen nucleus 



C = io- 



-i H - 



2.33m p corresponding to a helium 



abundance of 0.1 by number, and a magnetic field that fol- 
lows the standard tih 1 ^ 2 scaling for densities below 10 6 cm- 3 
(Myers & Goodman 1988) with a weaker dependence at 
higher densities: 



B 

ithG 



(n H /10 6 cm- 3 ) 1/2 ifn H < 10 6 cm- 3 , 
(rtH/10 6 cm- 3 ) 1 ' 4 otherwise. 



(22) 



When integrating over the MRN size distribution we have 
assumed that the ratio of charged to neutral grains is inde- 
pendent of grain radius. 

The upper panel of Figure |l| shows the charged parti- 
cle densities for the single-size grain model. Ions and elec- 
trons are the dominant charged species for tih ^ 10 10 cm" 3 , 
negatively-charged grains replace the electrons for tih ii 
10 11 cm -3 , and above 10 14 cm -3 ions are dominated by pos- 
itively charged grains. The three components of the conduc- 
tivity tensor are plotted in the lower panel, cry dominates 
the other two components until nn J> 10 15 cm -3 . Below 
about 10 10 cm" 3 , the Hall term oi is negative and is typi- 
cally an order of magnitude below the Pedersen term, a 2 ■ ci 
changes sign at roughly 10 10 cm- 3 . Above this density, \a\\ 
approaches 02 and the ambipolar diffusion approximation 
breaks down for nu 10 13 cm -3 . 

The contributions of different charged species to the 
components of a are plotted in Figure ^. It should be noted 
that each species contributes positively to both cry and 02, 
but that the contribution to a\ carries the sign of Zj. The 
conductivity parallel to the magnetic field is dominated by 
the electrons, even when much less abundant than the other 
charged species, because of their low rate coefficient for neu- 
tral scattering. As one might expect, the grain contributions 
are negligible because of their large collision cross-section 
with the neutrals. Negatively charged grains dominate the 
Hall term below 10 10 cm -3 at which point the ion Hall pa- 
rameter has dropped to order unity, the ion contribution 




log n H (cm" ) 

Figure 3. As for Fig. hi but for an MRN grain-size distribution. 



takes over and a\ becomes positive. At the highest densities, 
positive grains are also important. The Pedersen conductiv- 
ity 02 is determined by the ions, except at the very highest 
densities when electrons start to contribute. 

The abundances and conductivity for the MRN model 
are plotted in Figure [| There are many more grains than 
in the single-size grain model, and grains therefore play a 
significant role in the ionisation balance at lower densities 
(Nishi et al 1991). Negatively charged grains dominate the 
electrons for ?ih <i 10 7 cm -3 , and positive grains dominate 
ions above 10 10 cm- 3 . The drop in electron no. density re- 
duces cr || considerably, but cry is still much larger than the 
other components for tih Ss 10 12 cm -3 . a\ becomes compara- 
ble to o"2 between 10 7 cm -3 and 10 11 cm -3 . It drops sharply 
at higher densities, changing sign at 10 12 cm -3 . 

Fig. ^ shows which species contribute to the compo- 
nents of er. Electrons again dominate ay, but a\ and 02 are 
determined by the grains. Note that for tih ii 10 10 cm~ 3 , the 
positive and negative grains contribute nearly-equal terms 
of opposite sign to the Hall conductivity. The slightly de- 
pressed abundance of positive grains (of order 0.1-1%) is 
sufficient for the net grain contribution to a\ to dominate 
<T2. Above 10 11 cm" 3 the asymmetry drops to the point that 
<7i| is a small fraction of 02, becoming dominated by the 
ion contribution above 10 12 cm -3 . 

Finally, the effects of a PAH population are shown in 
Figures |H| and |^. Ions and PAHs are the dominant charged 



species between 10°-10 9 cm 3 



For n H £ lO^m" 3 PAHs 



and grains of either sign are important. Despite the change 
to the charged particle abundances, the components of the 
conductivity tensor behave similarly to the MRN-only case, 



with |<ti| « 02 between 10 and 10 cm 



Electrons con- 
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Figure 4. As for Fig. ^| but for an MRN grain-size distribution. 

tinue to dominate an below 10 10 cm -3 , with PAHs dominat- 
ing at higher densities. Grains determine <xi and 02, with 
PAHs becoming important above 10 10 cm -3 . 

A broad range of magnetic field strengths are appro- 
priate to any particular value of the gas density, so it is 
important to consider the effects of departures of the ma; 
netic field from the scaling (p^). Thus we show in Fig, 
how the ratio |o"i|/o"2 depends on both B and tih for the 
MRN grain model. Consider first the solid contours, which 
refer to a gas temperature of 30 K. The innermost contour 
corresponds to a ratio of unity, with successive contours in- 
dicating where |<Ti|/(T2 = 0.1, and 0.01. The gradients are 
weakest parallel to lines of constant B/nn as the Hall pa- 
rameters are constant and a varies only as a result of the 
density-dependence of the n^.The dashed trajectory indi- 
cates the scaling of B with iin from eq. (p^), thus Figure 
|§] shows a for cuts taken along this trajectory. It is clear 
that |<ti| 02 over a significant region at densities charac- 
teristic of dense cores. In addition there is a region at higher 
densities (tih ~ 10 11 -10 13 cm -3 ) and mG field strengths in 
which the the ratio is greater than 0.1. The dashed curves 
show the effect of increasing the temperature to 300 K, which 
increases the electron and grain momentum-transfer rate co- 




log n H (cm" ) 

Figure 5. As for Fig. [I| but for an MRN grain-size distribution 
and PAHs 



efficients and generally increases \<J\\j(Ji (note that the ef- 
fects of the increased temperature on the charged-particle 
abundances have not been included). 



4 ALFVEN WAVES 

The results of the previous section show that the Hall term 
is significant for molecular gas densities in the range 10 s - 
10 11 cm" 3 . As outlined in §j2| this will have consequences for 
dynamical processes, particularly for the time-evolution of 
magnetic fields from an initial state. By way of illustration, 
we consider here the propagation of Alfven waves in weakly- 
ionised media. 

Linear wave propagation in molecular clouds has pre- 
viously been studied by Kulsrud & Pearce (1969), who as- 
sumed a plasma comprised of neutrals, ions and electrons, 
and showed that short wavelengths are rapidly damped by 
ion- neutral friction. Pilipp et al (1987) examined the influ- 
ence of charged grains, assumed to be all of the same size. 
Pilipp et al considered interactions between charged species 
and the coupling between the charged and neutral grain flu- 
ids that results from charging and neutralisation processes. 
Propagation at an oblique angle to the magnetic field has 
been considered by Cramer & Vladimirov (1997). 

Here we present a simpler treatment using the formula- 
tion of which neglects the inertia of the charged species, 
and thus is valid for wave frequencies ui below the collision 
frequencies jjp of the charged species with the neutrals. This 
also permits us to neglect the coupling between charged and 
neutral grains, which generally is only important at even 
higher frequencies (Pilipp et al 1987). 
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We begin by estimating the frequency ui c at which 
flux freezing breaks down, and showing that lu c <C ypj- 
This guarantees that our treatment applies to the interest- 
ing regime uo ~ uo c . For long- wavelength disturbances, the 
cV X E' term in the induction equation (Q) is negligible and 
the effects of finite conductivity are small: the dispersion re- 
lation for Alfven waves is simply u ~ kvA- At successively 
shorter wavelengths the cV X E' term becomes increasingly 
important, and the wave modes are strongly modified by 
the finite conductivity when w ~ (lo/va) 2 c 2 /4~k<j±, i.e. at a 
critical frequency 



(23) 



per 



Substituting for a± using equations ([jj]) and (|l^), and not- 
ing that the norm of a sum of vectors is less than the sum 
of the norms, we obtain a limit on lo c : 



E 



E 



7jPj ■ 



(24) 



Strict equality is approached in the limit \/3j \ — > oo, i.e. when 
all of the charged species are tied to the magnetic field lines, 
in which case ui c is the collision frequency of the neutrals 
with any of the charged species. More generally, uj c is less 
than this and, as pj -C p, is clearly much less than 'ypj 
for any charged species j. The smallest jjp is that for the 
largest (i.e. 0.1 pm) grains. In dense clouds, jgPg ~ JiPi, 
thus uj c ~ "iapG- As the mass in grains is of order a per- 
cent of the neutral gas, and only a fraction of the grains are 
charged, this implies that the dispersion relation we obtain 
below is valid for frequencies up to at least 100 lj c . At higher 
frequencies the neutral component does not play a signifi- 
cant dynamical role but merely provides a background drag 
that damps small-scale disturbances in the ionised compo- 
nent of the fluid. 



Figure 6. As for Fig. g but for an MRN grain-size distribution 
and PAHs. 
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Figure 7. The ratio \o\\/<?2 for temperatures of 30 and 300 K 
(solid and dotted contours respectively). The contours are spaced 
logarithmically, with levels 1, 0.1, and 0.01. The dashed curve 
indicates the magnetic field scaling adopted in calculating Figs |l] 



4.1 Linearisation 

We linearise equations @, (g), (§), (|), and @ about a 
homogeneous undisturbed state with v — 0, J = E' = 0, 
and B — Bx. Changes in the conductivity associated with 
the perturbations do not appear in the linearised equations, 
as E' vanishes in the unperturbed state. Thus we need not 
make any assumptions about how the charged particle abun- 
dances respond to the perturbations - all that is needed is 
the conductivity tensor of the unperturbed fluid. 

Considering waves propagating parallel to B of the form 
exp i(uit — kx) , and discarding an additional quadratic factor 
that yields the familiar isothermal sound waves propagating 
parallel and antiparallel to the magnetic field, we find that 
uj and k are related by the two dispersion relations 



lo 2 =p k 2 v\ exp(±i6 



j/u 



, 2 2 

k v A 







where 



_1 (kil/ax) , 



(25) 



(26) 



may range between and n/2. The two dispersion relations 
yield four roots; assuming that k > we are interested in 
the two roots with positive real part. Setting 



b = ^— ^- exp(H0) 
2u) c 



(27) 
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these two roots may be written 
wi = (y/T+lP + b) kv A 
and 



LU2 



i 2 2 

k v A 



(28) 



(29) 



where is the complex conjugate of wi. Note that |oia| < 
kvA < 

The conductivity parallel to the field, cry does not ap- 
pear in the dispersion relation because E'n = for these 
modes. The conductivity perpendicular to B, a±, deter- 
mines ui c and sets the length scale at which the dispersion 
relation is modified from the ideal MHD case by the finite 
conductivity of the plasma, and the relative sizes of the Hall 
and Pedersen conductivities sets 8. 

Solving the linearised equations for the perturbations, 
we find: 



5B X = 0, SB y = ±i sign(<7i)(5B z 



SJ 



±sign(cri)^<5S . 



Sv _ (kv A \ SB 

VA \ U) ) B 

SE' = -isign(cr 1 )exp(±i6l) — — SB . 

UJ C c 

The perturbation to the electric field in the rest frame of the 
unperturbed fluid is 



(30) 
(31) 

(32) 
(33) 



SE = ^s ig a(a 1 )(^-) V fSB. 



(34) 



Equation ( |3C| ) shows that the wave modes are circularly po- 
larised, with the sense of polarisation depending on the sign 
of (J\. 



4.2 Wave propagation 

The wave frequency uo is plotted as a function of k in Fig. 
^| for different choices of |cti|/<72- The solid curves give the 
real and imaginary parts of u for o\ = 0, the pure ambipolar 
diffusion case (see Kulsrud & Pearce 1969). For kvA/^c < 2 
the two circularly-polarised modes have the same frequency, 
thus one can construct linearly-polarised Alfven waves. At 
low wave numbers, the modes reduce to Alfven waves in the 
combined neutral and ionised fluid. As the wave number is 
increased the modes propagate at slightly below the Alfven 
speed and are damped, the damping rate increasing as k 2 . 
The phase speed drops significantly below the Alfven speed 
at kvA ~ w c . At kvA/^c = 2 the waves become evanescent, 
with two distinct modes. 

The other curves plotted in Fig. |§] show how the modes 
are modified as the Hall conductivity is increased in impor- 
tance. In general the damping rate is decreased because the 
Hall current is perpendicular to E' and therefore does not 
contribute to the energy dissipation in the fluid. When the 
Hall conductivity is non-zero, the evanescent regime disap- 
pears. 

The two circularly-polarised modes propagate with dif- 
ferent phase speeds and damping rates because the Hall con- 
ductivity introduces a handedness into the fluid - the micro- 
physical asymmetry between the positively- and negatively- 
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Figure 8. The dispersion relation for long-wavelength waves 
propagating parallel to the magnetic field in a uniform, weakly- 
ionised plasma. The real and imaginary parts of the wave fre- 
quency u) (upper and lower panels respectively) are plotted as a 
function of (real) wave number, k. The characteristic frequency 
uic is defined in eq. The curves are distinguished by the ratio 

of the Hall and Pedersen conductivities, | cri | /tT2 : solid — 0; dotted 
— 0.1; short-dashed - 1; long-dashed - 10. For clarity, the curves 
for the real part of ui for | o\ \ /era = 1 , and for the imaginary part 
for |ci|/(T2 = 0-1 have been omitted as they are similar to those 
for |<ti|/<T2 = 10 and respectively. 



charged species is manifested as a dependence on the sign 
of the magnetic field, or the sense of circular polarisation. 

For kvA/uc 3> 1, the roots of the dispersion relations 
with positive real part are 



and 

iOc 



kVA 



exp(i^) 



exp(i#) 



(35) 



(36) 



dB/dt 




LO 


V x (w x B) 




kVA 



A comparison of terms in the linearised induction equation 
yields 



(37) 



For u)\ or u>2 this ratio is (kvA/uJ c ) 2 or (kvA/^c)~ 2 respec- 
tively. Thus V X E' in the induction equation is balanced 
by dB/dt and V X (v X B) respectively. Equation ( p2"| ) 
shows that for these two modes, the ratio \SB / B\/\Sv / va\ is 
(kvA/oJc) or (fciu/^c) -1 ■ The highly-damped mode Wi mode 
can be regarded as a wave in the ionised component fluid 
that is damped by collisions with the nearly static neutrals. 
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The 0J2 mode corresponds to transverse oscillations in the 
neutrals with small perturbations in the magnetic field and 
ionised component (c.f. Kulsrud & Pearce 1969). Note that 
both modes are undamped in the limit that 02 — > 0. 

In the low- frequency limit kvA/^c < 1, Wi and ui2 be- 
come, to first order in kvA/w c , 



lo , . kv A n .kv A . , 
~ 1 ± — — cos u — % — — sm f 

kVA ^UJc 2uj c 



(38) 



Thus the real part of lu (i.e. phase speed) is hardly affected 
by a, but the damping rate is dependent on both ai and 02. 
It should be noted that including the Hall term decreases the 
damping rate for a given k by increasing uo c and introducing 
a factor of sin#. Thus the damping rate is reduced by a 
factor of sin 2 9 once <J\ is taken into account. 



5 DISCUSSION 

The evaluation of the conductivity tensor presented in §^ 
shows that for an MRN grain-size distribution the Hall 
conductivity is important for densities between 10 7 and 
10 11 cm -3 . PAHs, if present, reduce this range to 10 9 - 
10 10 cm -3 . The magnitude of o\ is sensitive to the slight 
difference (1% or less) in the number densities of positive 
and negative grains, and positive and negative PAHs, so a 
careful calculation of the ionisation balance at these den- 
sities is required to confirm these results. The sensitivity 
of the conductivity tensor to small changes in the relative 
abundances of these species potentially couples the qualita- 
tive dynamical behaviour of dense molecular gas to chem- 
ical reactions that affect the ionisation level in the gas, as 
these in turn may depend on the relative drift of the charged 
species and neutrals. 

In dense molecular clouds the grain-size distribution 
is modified by the agglomeration of small grains through 
thermal collisions and the sweeping up of smaller particles 
by large particles (Ossenkopf 1993; Ossenkopf & Henning 
1994). These processes remove most of the grains that have 
Hall parameters of order unity (grain sizes of order 100 A), 
which tends to reduce the size of |oi|. In this case one ex- 
pects the conductivity to resemble that for the 0.1 micron 
grain case (see Figs 1 and 2). 

At a given density and magnetic field strength, the con- 
ductivity also depends on the assumed cosmic-ray ionisation 
rate per unit volume, C n H- There is a simple scaling law for 
a as the fractional abundances nj/nn of the charged species 
depend only on the ratio jih/C (e.g. Nishi et al 1991). This 
allows a to be written as a function of Uh/C an d B/C rather 
than of tih, B and £ independently. Thus, for example, if 
£ is increased from our adopted value of 10 _17 H _1 s _1 to 
10 — 18 H _1 s _1 , the contours in Fig. (?] are shifted upwards 
and to the right by 1 logarithmic unit in each direction, sig- 
nificantly increasing the area covered by the |<7i| = 02 con- 
tour. The cosmic-ray flux is substantially reduced for densi- 



ties in excess of 10 13 cm -3 , which occur near the midplane in 
the inner few AU of protostellar disks, because the cosmic- 
ray flux is substantially reduced when the disk surface den- 
sity exceeds the cosmic-ray attenuation column, 96 gem -2 . 
In particular, the surface density at 1 AU from the Sun in 
the minimum solar nebula is approximately 1200 gem -2 
(Hayashi 1981). Under these conditions the charged particle 



abundances are substantially reduced and the Hall contri- 
bution is negligible. 

The Hall term introduces a definite handedness to the 
dynamical behaviour of the gas because of the differences in 
the properties of negatively- and positively-charged species. 
There is little change in dynamics on the largest scales, for 
which flux-freezing is a good approximation. On small length 
scales, the dynamical behaviour of the fluid is qualitatively 
altered, and is no longer unaffected under a global reversal 
of magnetic field direction. This has been previously noted 
in models for the acceleration of outflows from protostellar 
disks (Wardle & Konigl 1993), where the relative orientation 
(i.e. parallel or anti-parallel) of the large-scale threading the 
disk and the angular velocity of the disk affects the ability 
of the field to centrifugally drive a disk wind. These effects 
are also apparent in models of C-type shock waves (Wardle 
1997). 

We chose to illustrate these effects by examining the 
consequences for Alfven wave propagation, as the wave 
modes encapsulate the dynamical response of the fluid to 
perturbations. On large length scales, the Alfven modes are 
damped more rapidly if the Hall contribution to cr is in- 
cluded, which could in principle have consequences for the 
survival of the Alfven waves believed to support molecular 
clouds parallel to the magnetic field (Arons & Max 1975; 
Fatuzzo & Adams 1993; Gammie & Ostriker 1996), and 
for the contribution of wave dissipation to cloud heating 
(Zweibel & Josafatsson 1983). However at the relevant den- 
sities |<ti| <C (72, and in any case the factor of two differ- 
ence that would obtain when \cri\ ~ 02 is of the same or- 
der as the theoretical and observational uncertainties in the 
field strength and charged particle abundances that deter- 
mine the damping rate. Of more importance are the changes 
that occur on the length scales on which flux-freezing breaks 
down, when the two circularly-polarised Alfven modes have 
different phase speeds and damping rates. For example, con- 
sider the following tho ught experiment that follows from the 
linear analysis of §4.1. Imagine at some initial time a homo- 



geneous, static fluid permeated by a magnetic field of the 
form B — Bqx + B\ sin kxy, where So and Bi <C Bo are 
positive constants. How will the system evolve with time? 
In the ambipolar diffusion limit B y will oscillate back and 
forth with nodes where sin kx = 0, the oscillatory motion 
being damped by ambipolar diffusion. If the Hall term is 
important the motion is more complicated as one of the cir- 
cular polarisations is more rapidly damped than the other. 
The field tends to twist back and forth, with an associated 
sloshing of the fluid. Thus the evolution is quite different in 
the two cases. 

As the transport of angular momentum by magnetic 
fields and the breakdown of flux-freezing are believed to 
play important roles in the collapse of cloud cores to form 
stars (e.g. Shu, Adams & Lizano 1987; Mouschovias 1987), 
and the densities pass through the regime in which the Hall 
term is important, the consequences for the dynamics of col- 
lapsing cloud cores may be profound. Norman & Heyvaerts 
(1985) also noted that the Hall component of the conduc- 
tivity tensor may become significant during star formation. 
However they dismiss its effects, claiming that the Hall cur- 
rent would rapidly lead to a slight charge separation, with 
the associated electric field quenching the Hall current. Al- 
though this can certainly occur for particular idealised ge- 
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ometries with well-defined boundaries, it is not clear that 
this mitigation of the Hall effect is generally true. For ex- 
ample, consider the quasistatic collapse of an axisymmetric 
cloud core through a magnetic field, and suppose that at 
some initial time the core is not rotating and the magnetic 
field is poloidal. Then the current is toroidal and the electric 
field in the neutral fluid frame is poloidal. Thus X7xE, hence 
dB/dt, is poloidal, as is the Lorentz force J X B/c. Thus 
the fluid velocities and the magnetic field remain poloidal 
as the system evolves. In contrast, if the Hall term is im- 
portant, at t = the electric field has a poloidal component 
(as J is still toroidal initially) and will grow with time, 
forcing the fluid to also pick up a toroidal component. It is 
not hard to show that if the Hall and ambipolar diffusion 
components of the conductivity are of the same order that 
the toroidal components of the velocity and magnetic field 
will contribute significantly to the gravitational support of 
the core and therefore will significantly affect the dynam- 
ics of gravitational collapse. In this example the boundaries 
on which the charge is presumed to build up initially would 
have to be surfaces of constant azimuth, which cannot exist 
under the assumption of axisymmetry. 



6 SUMMARY 

In this paper we evaluated the conductivity tensor for 
weakly-ionised, dense molecular gas assuming a temperature 
of 30 K and a cosmic-ray ionisation rate of 10 _17 s _1 H _1 . 
We also considered the propagation of Alfven waves in this 
medium. Our results are summarised as follows: 

(i) For single-size grains of size 0.1 /xm, the Hall term 
does not become significant. 

(ii) For an MRN grain-size distribution the Hall term be- 
comes important between 10 7 and 10 11 cm -3 . 

(iii) PAHs, if present, reduce this range to 10 9 -10 10 cm . 

(iv) The magnitude of the Hall conductivity is sensitive to 
the small differences in abundances of positive and negative 
grains, and positive and negative PAHs. 

(v) Long-wavelength Alfven waves, for which flux freezing 
almost holds, are not much affected by the Hall term, apart 
from a reduced damping rate. Generally this reduction is 
not significant given the uncertainties in the parameters of 
the gas in molecular clouds. 

(vi) At wavelengths comparable to or less than the length 
scale on which flux-freezing breaks down, left- and right- 
circularly polarised Alfven waves of the same frequency have 
distinct phase speeds and damping rates. For the typical 
case, the right-hand (left-hand) polarisation is most strongly 
damped for propagation parallel (anti-parallel) to the mag- 
netic field. The other mode is damped on the ion-neutral 
coupling time scale, even at the shortest wavelengths. 

We noted that the modification to Alfven wave propagation 
is indicative of the change in the dynamical behaviour of 
the gas. In particular we suggested that the dynamics of the 
collapse of dense cloud cores may be profoundly affected. 

We thank Leon Mestel for useful discussions. C. Ng was 
supported by an SRCfTA vacation scholarship. The Special 
Research Centre for Theoretical Astrophysics is funded by 
the Australian Research Council under the Special Research 
Centres programme. 
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